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Abstract 

The knotting probability is defined by the probabihty with which an A^-step self-avoiding 
polygon (SAP) with a fixed type of knot appears in the configuration space. We evaluate these 
probabilities for some knot types on a simple cubic lattice. For the trivial knot, we find that the 
knotting probability decays much slower for the SAP on the cubic lattice than for continuum 



models of the SAP as a function of N. In particular the characteristic length of the trivial knot 
that corresponds to a 'half-life' of the knotting probability is estimated to be 2.5 x 10^ on the 
cubic lattice. 



1 Introduction 

The Self- Avoiding Polygon (SAP) with fixed topology gives a simplified model of real ring 
polymers in solution that have a topological constraint as well as excluded volume. Throughout 
the time evolution, a circular polymer keeps the same knot which is given to it when it is made; 
it does not change its topology under any thermal fluctuations since no crossing through itself 
is allowed. On the other hand, SAP corresponds to the special case of the Self-Avoiding Walk 
(SAW) that returns to the origin. If we construct a set of SAPs, then their topological states 
may contain several different knots. Therefore, it is not trivial how to realize the topological 
constraint on an SAP. One possible method for assigning the topological constraint on an SAP 
is that after generating a large number of SAPs we select only such SAPs that have the same 
given knot. By this method, we can assign the topological constraint on any model of SAP. In 
this context, the probability that a given SAP has a fixed knot plays a central role, and we call 
it the knotting probability of the model of SAP for the knot. Among many different models, the 
SAP on the cubic lattice with fixed topology is one of the most fundamental models of SAP. It 
has an advantage that the definition is very simple. We expect that the model should be suitable 
for general and mathematical study. In fact, several rigorous results on knotting probability 
have been derived for the SAP on the cubic lattice [||, ^ Thus, the main motivation of 

the present research is to characterize the knotting probability of the SAP on the cubic lattice 
through numerical simulations. 

Let us discuss some previous numerical results on knotting probabilities of SAP ^, |7|, 
|, |, 0, |n|, |13|, For a model of SAP (or random polygon) with N steps, we 



denote by Pk{N) the knotting probability of the model for knot K. For the Gaussian model of 
random polygon and the rod-bead model of SAP, knotting probabilities were evaluated through 



numerical simulations for the trivial knot K = % ^ and also for some non-trivial knots 

|l^, |l3|, |l^, 15 1 . It was found that for the models of SAP and random polygon, the knotting 
probability as a function of the step number N is given by the following: 

Pk{N) = C{K)N'^^^^ eM-N/N{K)l (1.1) 



where C{K), m{K) and N{K) are fitting parameters. For large A^, the formula is consistent 



with the asymptotic expansion of the partition function of the SAP with the fixed knot type. 



For finite A^, although it is not rigorous, it seems that the formula (1.1) fits numerical data 
well. From the numerical result it was conjectured that the parameter N{K), which we call 
the characteristic length of knot K, should be given by the same value for any knot K [^]. 
Furthermore, it was also conjectured that the parameter m{K), which we call the exponent of 



knot K, should be universal for different models of SAP or random polygon 14, 15 1. We 



note that the fitting formula (1.1) together with the two conjectures are consistent with the 
standard asymptotic behavior expected for SAP or random polygon. Here we also note that 
the rod-bead model is an off-lattice model of SAP. 

Recently, the SAP on the cubic lattice with fixed knot was studied through a numerical 
simulation using the BFACF algorithm, which generates SAPs with the same fixed knot but 
different step numbers N |]l^]. In the simulation, the exponent m{K) and the growth constant 
for the number of allowed configurations of the SAP with knot K has been estimated for some 
knots. Furthermore, it was shown that the knotting probability for the trivial knot decays 
"exponentially" for a face-centered cubic lattice |11] and for the cubic lattice However, any 



precise estimate of the knotting probability or the characteristic length N{K) has not been given 
for the cubic lattice. Thus, it is the primary purpose of this paper to evaluate the characteristic 
length N{K) for the SAP on the cubic lattice. 

In 1962, Delbriick [|l^ noticed that the topological constraint may be very important for 
polymers in biology and chemistry. Since then, the topological problem has been studied in 
several approaches in physics and mathematics. As one of the earliest studies, des Cloizeaux 



and Mehta estimated through numerical simulations the critical exponent v for the internal 
distance of the Gaussian random polygon and discussed some possible properties of random 
polygons under the topological constraint. 

After the rediscovery of the pivot algorithm, many properties of SAW and SAP have been 
investigated not only in the field theory but also by computer simulations |2l|. It seems, 
however, that there are only a few works such as 19| where the gyration radius {I^)sAP of 
the SAP on the cubic lattice is studied by numerical simulations. Let us denote by {R^)sAW and 
{R'q)sap the mean square of the gyration radius of SAW and SAP, respectively. It is interesting 
to evaluate the universal amplitude ratio between {Rq)sap and {Rq)saw- The ratio has been 
evaluated only up to 0{e) in the e-expansion method [^]. Furthermore, according to the scaling 
theory of polymers, the exponent z^sap for {Rq)sap should be given by the exponent t'sAW for 
{Rq)saw- The agreement is confirmed up to 0(e) by the renormalization group theory p^j . 
However, it is nontrivial to confirm the agreement for the SAP on the cubic lattice through 
numerical simulations. Thus, the numerical study of the gyration radius for the SAP on the 
cubic lattice is another purpose of this research. 

Hereafter in this introduction, we explain some of the main results of our numerical sim- 
ulations. Employing the pivot algorithm, we construct a large number of SAPs for the SAP 
on the cubic lattice with given step number N. For the gyration radius, we have obtained the 
exponent: t'sAP = 0.5867 =b 0.0017. This is indeed in good agreement with the estimate of the 
critical exponent of the SAW in the e-expansion: z^sAW = 0.5882 it 0.0011. Thus, our simulation 
in the paper also confirms the agreement of the two exponents. 

Let us explicitly consider the method for the topological constraint on the SAP in the cubic 
lattice with given step number N. The pivot algorithm of the SAP can generate all the allowed 
configurations of the SAP with equal probability. Therefore, the set of SAPs generated by 
the algorithm may contain various knots. Suppose that we have constructed M SAPs of step 
number N. Calculating some knot invariants for each of the SAPs, we effectively detect the knot 
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types of the SAPs. We enumerate the number of such SAPs that have the same set of values 
of the knot invariants for knot K, and denote it by M{K). The expectation value of a physical 
quantity under the topological constraint with fixed knot K can be effectively calculated by 
taking the statistical average of the quantity for the M{K) SAPs. 

We now turn to the knotting probability. Let the symbol Pk{N) denote the knotting 
probability of the SAP with steps on the cubic lattice for knot K. If there are M{K) SAPs 
in the total M SAPs, then we evaluate it by Pk{N)= M{K)/M. In our simulation, 10^ SAPs 
(M = 10^) are constructed for six different values of the step number from N = 500 to = 3000. 
We have found that almost all SAPs are topologically equivalent to the trivial knot, and also 
that the resulting values of the knotting probability for the trivial knot are fitted well by the 
formula ( |1.1D . Thus, for the trivial knot, we have obtained an estimate of the characteristic 
length 

A^(0) = (2.5 ± 0.3) X 10^ (1.2) 

This result means that trivial knots are dominant among SAPs on the cubic lattice when the 
step number is less than 10^. It implies that when A^ > A^(0), the majority of SAPs on the 
cubic lattice have some non-trivial knots. The large value of the characteristic length might be 
a consequence of strong self-avoiding effect of the SAP on the cubic lattice. 

This paper is organized as follows. In section 2 we will estimate the universal amplitude ratio 
and the critical exponent z^sap of the gyration radius. In section 3 we will explicitly discuss the 
knotting probability for the trivial knot, and obtain the estimate of the characteristic length. 
We will also show that the knotting probability for the trivial knot decays almost linearly since 
the characteristic length of the trivial knot is so large. Finally, we will discuss a possibility that 
the characteristic lengths take an unique value without depending on knot types, in the last 
section. 
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2 The mean-square of the gyration radius 
2.1 Previous results on SAW and SAP 

An A^-step Self- Avoiding Walk (SAW) w in is a sequence wo,wi, - ■ ■ ,wn of N + 1 distinct 
points in such that each point Wi is one of the nearest neighbors of its predecessor Wi^i: 
{wi — Wi-il = 1 for i = 1, • • • , N. It is also subject to a constraint that any site never be occupied 
by two or more points in a sequence {wi}. The points wq and wn are the endpoints of w. The 
components of Wi are represented by w^"^ for a = 1,2,3. The Self- Avoiding Polygon (SAP) is 
a special case of the SAW that makes a ring. We consider an (N — l)-step SAW and denote it 
hy w = {wq^wi, ■ ■ ■ jWn-i}- If the endpoint wq is the nearest neighbor of the endpoint wn-i, 
this is an A^-step SAP in Z^. 

The mean-square of the gyration radius of the SAP {Rq)sap is smaller than that of the 
SAW {Rq)saw- This fact comes from the following difference. The endpoints of SAW are free, 
while those of SAP are constrained, that is, the endpoints of SAP should meet at the nearest 
neighbor sites of the lattice. This means that the SAW has a possibility of a longest end-to-end 
distance whereas the SAP does not. 

The RG argument gives several results of the SAW and the SAP. Among these, we are 
most interested in the amplitude ratio {R'q)sap/{Rg)saw- According to the RG theory, the 
amplitude ratio should be universal, i.e. it does not depend on the details of the models for 
the SAW or the SAP. The ratio has been evaluated by using the RG equation and the cluster 
expansion up to 0{e) 23 1: 

{RDSAP 



0.568. 2.1 



GlSAW 



The ratio has also been estimated in Ref. |19] through numerical simulations with N < 800: It 
is given by 0.538 ± 0.006. 

Another interesting result of the RG argument is that the critical exponent for the mean- 
square of the gyration radius is universal. It is well-known that the exponent fsmi of SAW 
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corresponds to the critical exponent of the 0(n) vector model in the limit of n going to zero. 
The precise estimate of the critical exponent z^saw has been made by the e-expansion method 



through this correspondence [EJ, QB, 26, 



zyg^^ = 0.5882 ± 0.0011. (2.2) 

The exponent i^saw has also been precisely evaluated by the Monte Carlo simulations, where 
the best estimate is given by the following ||2C 



= 0.5877 ± 0.0006. (2.3) 

2.2 The amplitude ratio 

In our simulation of SAPs on the cubic lattice, we have employed a length-conserving dy- 
namical algorithm which keeps endpoints fixed. This algorithm was investigated in detail by 



Madras, Orlitsky, and Shepp |20, 21|. We call this algorithm MOS pivot for short. Some details 
are described in Appendix A. 

We calculate the ratio {Rq) sap / {Rg) saw and the exponent i^sap using the MOS pivot and 
compare resulting values with theoretical and other simulated values. The mean-square of the 
gyration radius is calculated by 

T 



{Rg)sap = Tf;^ 



T 

t=i 




(2.4) 



where T is the number of polygons and the step number, and the symbol ivf^ denotes the i-th 
site in the t-th SAP in the sequence of SAPs given by "200 successful MOS pivot operations" . 
(See Appendix A for details.) 

Let us note that the gyration radius should have the asymptotic behavior 

{RDsap = AsApN^"^^^ (1 + BN-^) as iV ^ oo. (2.5) 

Choosing ^ for the exponent A of the correction term, the numerical data given by the equation 
(O) are plotted in figure | and fitted to ^sapA^^''^'"'(1 + BN'^), where ^sap, B and i^sap are 



fitting parameters. Since there are not enough data to fit in a four-parameter curve in our 
simulation, we have assumed the fixed value for the exponent A. We plot the mean-square 
of the gyration radius from N =500 to N =30000 in figure I Here, we set T, the number 
of polygons, to 10000. The error bars denote one standard deviation given by the Poisson 
distribution to the number T of polygons with N. 

The three parameters are obtained by fitting the data in figure ^ 

z^sAP = 0.5867 ±0.0017, (2.6) 
^SAP = 0.1101 ±0.0037, (2.7) 
B = -0.06 ±0.04. (2.8) 

The errors are subjective 68.3% confidence intervals. For the SAW, we recall that Madras et 
al. [^] estimated the exponent z^saw by the pivot algorithm up to N =80000. The estimate of 
J^SAW is given by ( p. 3]) together with the following estimates of the fitting parameters 

^SAW = 0.19455 ±0.00007, (2.9) 
B' = -0.11432 ± 0.00465, (2.10) 
A = 0.56 ±0.03. (2.11) 



Here, B' corresponds to AB in our estimated values. The estimate ( p.6[ ) is in good agreement 
with that of the RG argument ( ^.2| ) and that of the simulation for the SAW (|2.3| ). 



From (2^) and ( |2.9| ), we obtain the amplitude ratio 

j^J^^ = 0.566 ± 0.019, (2.12) 

{Rgisaw 

where an error is a subjective 68.3% confidence interval. The amplitude ratio ( p.l2| ) is consistent 
with the estimated value of the RG argument ( p.l]) . The estimate of Ref . ||l^ is a littile smaller 
than (2.12). However, it should be consistent with the estimate of Ref. |19|, if we consider the 



fact that the step numbers N of the simulations in Ref. [|19[ are much less than 30000. 
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3 The characteristic length A^(0) 

3.1 A method of evaluation of the knotting probabiUty 

We evaluate the knotting probability in the following way. We generate M SAPs of N steps 
and then enumerate the number M{K) of those polygons which are equivalent to a given knot 
type K. We define the knotting probability Pk{N) by the ratio M{K)/M. 

In our simulation, we determine knot types of SAPs using the 2nd-order Vassiliev-type 
invariant and the Alexander polynomial evaluated at t = —1. The Vassiliev-type invariants 
have the following advantages: (1) we can calculate them in polynomial time, and (2) we can 
calculate them without consuming a large memory area |^^. The Vassiliev-type invariants are 
not complete invariants. However, in the practical sense we can safely say that if the value of 
the Vassiliev-type invariant computed for an SAP is zero, this SAP is a trivial knot. We will 
see in section |3.3| that complicated knots are very rare events in our data. Even if non-trivial 
knots are misidentified as the trivial knot by the Vassiliev-type invariant (this chance is very 
small), it would not affect the results of this paper. 

For calculating the knotting probability, we generate random sequences of SAPs. We make 
five seeds for each of the six step numbers N =500, 1000, 1500, 2000, 2500, and 3000, and then 
apply an operator P on them for a large number of times (20000 times). Here, P denotes 200 
successful MOS pivots. (See Appendix A.) Then, we get effectively random sequences of SAPs. 
A sequence of SAPs derived from each seed has a set of 20000 SAPs. Here the sequence consists 
of 20000 X 200 successful MOS pivot operations. For each step number A^, we thus get samples 
of 100000 SAPs. 

In order to analyze the behavior of the knotting probability, we use the fitting formula ( |1.1| ) . 
Here, we write it again: 

Pk[N) = C{K)N'^^^^ exp 
where C{N), N(K), and m{K) are fitting parameters. In particular N{K) is called characteris- 
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N{K) 



tic length with knot type K. This formula was introduced by Deguchi and Tsurusaki |12, 14, ^ 



They pointed out that the formula (IT) is suitable for the knotting probabilities of the Gaussian 
and rod-bead models. We will show that this is suitable also for the knotting probability Pk{N) 
of the cubic lattice model. 

3.2 The random events of knotted polygons 

Let us discuss statistics of knotted polygons generated in our simulation. We will see in sec- 



tion 3.3 that almost all polygons are of trivial knots. Therefore, we may assume that non-trivial 
knots are generated such as the Poisson random events: the number of trivial knots between 
two knotted SAPs will follow the Poisson distribution, if the SAPs are randomly constructed. 

Let us consider the "time interval" of the Poisson random events. We recall that the time 
t is a discrete number. We measure the length L of the time interval of SAPs from time ti 
at which a knotted SAP appears to time t2 at which the next knotted SAP appears after ti, 
and set L = t2 — ti. If non-trivial knots are generated as the Poisson random events, the time 
interval L follows the function: 

D{N, L) = A{N) P^{N)^-^. (3.1) 

We call it the discrete distribution function of the time interval L. Here, P^{N) is the knotting 
probability for the trivial knot. We have introduced A{N) for a technical reason and hence 
D{N,L) is not necessarily normalized. 

The discrete distribution functions of the time interval L are numerically evaluated as follows. 
The time interval L is a discrete random variable. We introduce a sequence of natural numbers 
{lo,li, ■ ■ ■ }, where /j = 50 x i. Then, we count the number of polygons with L taking values in 
[li, k+i) and plot it at /j for each i. This is the discrete distribution function of the time interval 
L obtained from numerical evaluation. 

We plot the discrete distribution functions obtained by the numerical evaluation and D(N, L) 
as a function of L, where P^i^N) is estimated from table |^ and ^(A^) is chosen to fit these 
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distributions (See figure 2a and 2b). Error bars denote one standard deviation. They are 
estimated by applying the Poisson distribution to the number of samples included in interval 
L. For N =1000, 1500, 2000, 2500, and 3000, these graphs show a fairly good agreement with 
the function D(N,L). In the case of N=500, the data deviate from the function. This is not 
unexpected since we have too few samples of knotted polygons. We do not plot the graph for 
N =500 in this paper. 



We expect from figure 2a and 2b that the MOS pivot makes a uniform ensemble for knots. 
Thus, we can calculate the knotting probability and the characteristic length using the MOS 
pivot. 

3.3 The number of unknotted polygons 

The table || gives the number of each knot type with respect to N. Here, errors correspond 
to 68.3% confidence intervals. They are estimated by applying the binomial distribution to the 
number M{K) of polygons for knot K. In the table |l| we explain the notations: 0,3i,4i denote 
the trivial knot, the trefoil knot and the figure eight knot, respectively. The other knot types 
are denoted by etc. 

The table tells clearly that almost all the generated SAPs are trivial knots. M(3i) is also 
much larger than M(4i). The other knots (etc.) are nearly equal to zero. 



Step Numbers 


M(0) 


M(3i) 


M(4i) 


etc. 


M 


500 


99849 ± 25 


147 ± 24 


3 ± 3 


1 ± 1 


100000 


1000 


99640 ± 38 


344 ± 37 


9 ± 6 


7 ± 5 


100000 


1500 


99430 ± 48 


541 ± 46 


24 ± 10 


5 ± 4 


100000 


2000 


99208 ± 56 


752 ± 55 


27 ± 10 


13 ± 7 


100000 


2500 


98965 ± 64 


985 ± 62 


38 ± 12 


12 ± 7 


100000 


3000 


98787 ± 69 


1157 ± 68 


40 ± 13 


16 ± 8 


100000 



Table 1: the Number of Generated Knots 
The estimated errors correspond to one standard deviation. 



Next we focus on the knotting probability for the trivial knot, and plot Pid{N) as a function 
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of N (figure |3|). The error bars are one standard deviation. P${N) decays linearly with respect 
to N. It is expected that -P0(A^) becomes decays exponentially when goes to infinity. In 
|]|, ^, ^ it was shown that the knotting probability P(ij{N) tends to zero "exponentially": 

P0(iV) = Aexp[-nN + o{N)], when iV ^ oo. (3.2) 

Thus, this is a natural situation. 

The asymptotic shape ( p.2| ) is realized for the trivial knot in our case when the fitting 
parameters of the formula ( |1.1[ ) take the following values: |m(0)| <^ 1 and C(0) ~ 1. In fact, 
using the least-squares estimation, we find: 

C(0) = 1.0035 ± 0.0035 (3.3) 

m(0) = (-4.7 ± 5.7) X 10"^ (3.4) 

iV(0) = (2.5 ± 0.3) X 10^ (3.5) 

= 0.748 (3.6) 

Prob{x^ > 0.748) = 0.862 (3.7) 

(errors are one standard deviation). Here, the value (i.e. sum of square of normalized 
deviations from the regression line) can serve as a criterion of good fit. It should be distributed 
as with n — 3 degrees of freedom, where n is the number of data points in the fit. Prob{x^ > 
0.748) is the probability that would exceed the observed value, in this case 86.2%. This 



implies that the formula (LI) is suitable. It is remarkable that the characteristic length of the 
trivial knot is much larger than the value expected from the rod-bead model. 

There have been a few simulation studies on several lattices. In |]ll| the knotting probability 
for the trivial knot was calculated on a face-centered cubic (FCC) lattice. It was shown that 
assuming the two-parameter fitting formula: P(i){N) = C{%)e~'^^^^^ , the parameters were given 
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by 

a(0) = (7.6 ±0.9) X 10"^ 

C(0) = 1.0011 ±0.003 

= 2.7 

Probix^ > 2.7) = 0.44, 

where errors were one standard deviation. Our interpretation of a(0) taking 7.6 x 10~^ is that 
the characteristic length is 1.3 x 10^. Comparing the cubic lattice with the FCC lattice, the 
characteristic length is larger on the cubic lattice than on the FCC lattice. In it was also 
shown that the exponent a;(0) in the above form was (5.7 ± 0.5) x 10^^ on the cubic lattice. 
This corresponds to the characteristic length of 1.8 x 10^. We can expect that this is consistent 
with our estimated value within the error bars. Here we note that the connection of Ref. [|l^ ] 
shall be discussed in section Q and also note that the simulation results of Refs. p9| , |30| ] seem to 
contain some information on the characteristic length of the SAP on the cubic lattice. However, 
we are unable to derive any appropriate estimate from them. 

Let us return to our data. We estimate not only the characteristic length A^(0) ~ 2.5 x 10^ 
but also the exponent of a correction term m(0) ~ 0. These parameters show that the knotting 
probability for the trivial knot serves a "pure exponential" decay on the cubic lattice. This is a 
new result for the cubic lattice model. Thus, the simulation in the paper improves that of Q. 

We have several interpretations of the characteristic length N ((}>). Since A^(0) is so large, we 
found that the number of knotted polygons are much smaller than that of unknotted polygons 
in our simulation. We expect that the knotting probability for the trivial knot decreases to 
about 30% at = N {(/>). Thus, non-trivial knots become majority of SAPs when A^ is larger 
than A(0). In addition to this, we will see in section ^ that A(0) ~ N{K) for any knot type K 
is due to the fact that A^(0) is large. 

In figure |^, we plot the knotting probability for the trefoil knot P^^ (A) as a function of A, 

13 



where error bars are one standard deviation. The data points almost he on a straight hne. Since 
the data of the trefoil knot are only six points in the step number, we cannot fit the knotting 



probability P3^{N) to formula (1.1). We leave the selection of fitting parameters for the trefoil 
knot or more complicated knots as future investigations. 

The knotting probability i-3^(A^) may have a finite-size effect. When we set m(3i) = 1 as 
expected from |l^, the formula ( |1.1| ) does not match with Ps-^^N) in our data. Deguchi and 
Tsurusaki argued that a finite-size effect appeared in the knotting probabilities in [14|. When 



we plot the straight line fitted to P3^(A^), the line intersects with the x-axis at a positive value. 
This is the finite-size effect. We introduce offset parameter Nini and replace by = N — Nini 
in the formula ( |1.1D . When we fix m(3i) = 1, our rough estimation gives as Nini ~ 140 for the 
trefoil knot. On the other hand, we roughly estimate Nini ~ for the trivial knot. We expect 
that such finite-size effect would also appear in the knotting probabilities for the SAPs with 
more complicated knot types. It could be confirmed by generating SAPs with much larger A^. 

In section p.2| we calculated the gyration radius including all possible knots, and then es- 
timated the universal amplitude ratio and the universal exponent i^. Without classifying knot 
types, however, we can effectively consider only trivial knots when A^ < A^(0). The univer- 
sal amplitude ratio and the universal exponent are evaluated effectively for the trivial knots, 
although 2% of the SAPs are non-trivial knots and we neglect their infiuence. 

4 A consequence of 

the large characteristic length A^(0) 

The trivial knot dominates among SAPs on the cubic lattice when A'^ is less than 2.5 x 10^. 
Our interpretation is that the cubic lattice has so strongly the excluded-volume effect that it 
almost prevents appearances of knotted polygons. 

We have a conjecture that the appearance of a complicated knot is a rare event on the cubic 
lattice. According to [0, 0, non-trivial knots occupy a large number of configurations of 
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SAPs for N > N{9). We believe that the above situation is reahzed also on the cubic lattice. If 
the formula (|l.lD is a suitable form of the knotting probability and if the trefoil and figure eight 
knots on the cubic lattice behave like these of continuum models, the ratio M(4i)/Af (3i) for each 
step number N should depend only on the ratio C(4i)/C(3i). We recah that M(3i) » M(4i) 
for N < 3000 from the table [l|. Then, at iV ~ 2.5 x 10^ almost ah SAPs are expected to be 3i 
knotted polygons on the cubic lattice unlike SAPs on continuum models. 

We see that our estimated value A^(0) is related to the growth constant fi from the viewpoint 
of Orlandini et al. [p^] . The asymptotic behavior of the number of A^-step polygons cat is given 
by 

= aA^°-V^ (1 + 6A^-^ + o(Ar-i)) , ^4_^) 

where a and A are critical exponents. For fixed knot type K, it is believed that the number 
cn{K) of polygons with knot K should have a similar form: 

cn{K) = a{K)N'^^^'^-^ ^i{Kf (l + h{K)N-^^^^ + o(A^"^)) . (4.2) 

Then, the knotting probability Pk{N) is given by cn{K) / cn ■, and this implies that the char- 
acteristic length A^(0) relates to the growth constants /i and /i(0). We can estimate the ratio 
from the value of N{%. 

JL. ^ e'/Nm ^ 1 + (4 ± 2) X 10-^ (4.3) 
/x(0) 

where an error is two standard deviations. 

Let us discuss the independence of the characteristic length N{K) with respect to knot type 
K. Orlandini et al. fl^ calculated the growth constants directly and showed the following 
equality: 

/i(0) = iu{K) = 4.6836 ± 0.0038 for any knot type K, (4.4) 
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where an error corresponds to 95% confidence interval. In addition Guttmann estimated the 
growth constant [^] 

/i = 4.68393 ± 0.00002 (4.5) 
using exact enumeration and series analysis (an error is one standard deviation). From (|4] 



and (4^), we obtain the ratio 



- ~ 1 + (7±8) X 10-^ (4.6) 



/x(0) /x(K) 



where an error is two standard deviations. From (4.3) and (4.6), we expect that the characteristic 
length is independent of the knot type: A^(0) ~ N{K) 2.5 x 10^ for any knot type K. 

Although the independence of ij-{K) with respect to knot type K is pointed out by Orlandini 



et al. [16], we can also confirm it more precisely through our simulation. While Orlandini et 
al. [jl^ calculated the growth constants, we have no direct calculation of them. However, we 
can predict that the differnce between /i(0) and ii{K) is very small when we use the following 
inequalities: liminf A^~^ log CAr(ii') > /i(0) and limsup A^^""^ logCAf(i^) < for any knot type AT, 
which were proven in 0, |3|. If both liminf A^~^ \ogCN{K) and limsup A^~^ logCAr(X) exist and 

N-*co N^oo 

take the same value fJ-{K), then we have /x(0) < ij,{K) < These inequalities and the estimate 
limit the ratio fj,{K) / to 



1 



<(4±2)xl0~^ (4.7) 



This is a strong bound. Thus, we can expect that /i(0) ~ f^iK) for any knot type K. 
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Appendix 

A The method for generating SAPs 

Let us discuss the method for making SAPs in the paper. We first construct a seed SAP 
w, and then derive a random sequence of SAPs {w^^^ , w^'^^ , ■ ■ ■ , w^^~^^ , w^^\ ■ ■ ■ }, where = 
P{w^^~^^) and w^'^^ = w{0). Here the operator P will be defined later. 

We construct a 'seed' SAP, by combining two SAWs which have the same endpoints, in the 
following way: first we make an A^/2-step SAW using the myopic self-avoiding walk (MS AW) 
algorithm |32|, where N is an even integer; secondly we perform the MOS pivot transformations 



with respect to k = 0,1 = N/2. [^]; finally we concatenate the endpoints of the new and 
original SAWs respectively (figure P), and get an SAP with the step number N if it has no 
self- intersections. 

Let us discuss how to construct a random sequence of SAPs in our simulations. If one of the 
MOS pivot transformations changes a given SAP into a different SAP, this operation is called 
a successful MOS pivot operation. We consider a sequence of successful MOS pivot operations. 
Then, we define an operator P by 200 successful MOS pivot operations in the sequence. We note 
that the SAP obtained by a single successful MOS pivot is not independent from the original 
one: they are correlated. However, the correlation decays almost completely after 200 successful 
MOS pivot operations, which shall be shown in Appendix B. Thus, we may consider that for 
any given SAP w, P{w) is independent from w. 

For generating SAPs, we use the Mersenne Twister that is a pseudo random number genera- 
tor |3^. This algorithm has the following properties: (1) we can get many samples because the 
period is 2^^^'^^ — 1, (2) we treat high dimensional space (max 623 dimensions), (3) pseudo ran- 
dom numbers are generated fast and (4) we can use the memory efficiently. Thus, the Mersenne 
Twister is a high performance generator. 
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B The decay of correlations among the SAPs ob- 
tained by the MOS pivot operations 

In order to check the vahdity of the random sequence of SAPs constructed in the paper, we 
show exphcitly how the correlation among the SAPs decays after applying a number of MOS 
pivot operations. 

Let us regard the number r of successful MOS pivot operations as the time of evolution of 
the shape of SAP under the sequential MOS pivot operations: the seed SAP w a random 
sequence corresponds to w^O) for r = 0; w{t) is defined by p^iw) for r > 0. Here p have denoted 
a successful MOS pivot operation. Let us now define the correlation function for the structure 
of SAP with the step number in the following 



C{t) — J— 



ELi Eto^ [ (0) - ^ (0) 

where wf^^ (r) is the a-component of the i-th. site of SAP (or SAW) after r pivot operations and 
(•) denotes the statistical average. 

In figure ^ and we plot the correlation C(t) versus the time r in the case of = 
500 and 1000, respectively. Error bars show one standard deviation which are estimated by 
assuming that the data follows the Poisson distribution. In our simulation, we evaluated the 
correlations of SAPs with the step number N at 500 and 1000, and generated 10000 sequences 
(starting from 10000 different seeds) to take the statistical average. 

The decay rate of the correlation function C(t) is slower in the case of the MOS pivot 



than in the case of normal pivot algorithm (figure ^ and 6b). This is due to the difference 
in the numbers of independent pivot transformations. Note that the number of independent 
transformations for the normal pivot algorithm ^ is (i!2°'— 1 (the normal pivot algorithm is 
an algorithm with one endpoint free while the other fixed) , which correspond to the number of 
all the elements of d-dimensional orthogonal group on a hypercubic lattice, while for the MOS 
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pivot 2d{d — 1) + 1 in the d-dimensional hypercubic lattice. In the cubic lattice, the MOS pivot 
has 13 pivot transformations while the normal pivot algorithm has 47. 
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Figure 1: Mean-Square Radius of Gyration (-Rg) 
We give the log-log plotted graph of the mean-square of gyration radius versus the step number N . 
Here error bars denote one standard deviation. 
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Figure 2a: Discrete Distribution Function of Interval L of Non-Trivial Knots for N — 1000 

The number of occurring a knotted SAP with respect to the time interval L is represented by the 
Poisson distribution: D{N,L) = A{N)Pq{N)^ . Pq{N) is the knotting probability of the trivial knot 
with the step number N and A{N) is a normalization factor. We call it the discrete distribution 
function of the time interval L. In our simulation, we count the number of polygons with L taking 
values in [Zj, Zj+i) and plot it at Zj for each subscript i. Here, the sequence {Iq, Zi, • • • } is defined by 
li — 50 X i. The data are fitted to the distribution D(N,L). We measure the length of the time 
interval of SAPs from time ti at which a knotted SAP appears, to time t2 at which the next knotted 
SAP appears after ti, and denote it by L (= ^2 — ^i)- 
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Figure 2b: Discrete Distribution Function of Interval L of Non- Trivial Knots for = 1500,2000,2500 
and 3000 

We plot the distributions D{N, L) and the data of L taking values in [/j, /j+i) for the step number 
from N = 1,500 to = 3,000, respectively. The data are fitted to D{N, L) well. Therefore, the 
MOS pivot generates SAPs without biasing statistics of knots. 
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Figure 3: Knotting Probability Pi/,{N) 
We give the graph of the knotting probabiUty P(i,{N) for the trivial knot versus the step number N. 
Here error bars denote one standard deviation. The data behaves as a liner decay with respect to 
N. 
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Figure 4: Knotting Probability Ps^iN) 
We give the graph of the knotting probabihty P^^^N) of trefoils versus the step number A^. Here 
error bars denote one standard deviation. If we assume a finite-size effect for the SAPs of trefoils, 
we fit the data to the straight line which intersect x-axis at positive value. 
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Figure 5: Generated Seed SAP with Step Number 
An original SAW w generated by the MSAW algorithm with the step number N/2 transforms into 
another SAW using the MOS pivot with respect to A; = and I — N/2. We concatenate the edges of 
new and original SAWs, respectively. We get an A"-step SAP. 
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Figure 6a: Pivot Correlation C{r) at N = 500 
We plot the pivot correlations C(r) versus the time r. The MOS pivot and the normal pivot 
operations in the case of N — 500, respectively. Here, the r is the number of trials, and error bars 
correspond to one standard deviation. The correlation of the MOS pivot decays slower than that of 
the normal pivot. 
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Figure 6b: Pivot Correlation C(r) at = 1000 
Similarly to the Fig.^ we plot the pivot correlations C(r) versus the time r. The MOS pivot and 
the normal pivot operations in the case of = 1000, respectively. Here, the r is the number of trials, 
and error bars correspond to one standard deviation. At r = 200, we consider correlations C(r) as 
zero focusing on the MOS pivot correlations for A^ = 500 and 1000. 
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